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Abstract: We study the electromagnetic on-shell form factor of quarks in massless perturbative 
QCD. We derive the complete pole part in dimensional regularization at three loops, and extend 
the resummation of the form factor to the next- to-next- to-leading contributions. These results are 
employed to evaluate the infrared finite absolute ratio of the time-like and space-Uke form factors 
up to the fourth order in the strong coupling constant. Besides for the pole structure of higher-loop 
QCD amplitudes, our new contributions to the form factor are also relevant for the high-energy 
limit of massive gauge theories like QED. The highest-transcendentality component of our results 
confirms a result recently obtained in 0\C=4 Super- Yang-Mills theory. 
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1. Introduction 



The electromagnetic form factor of quarks is a quantity of considerable interest in Quantum Chro- 
modynamics (QCD) and in gauge theories in general. At high photon virtualities this quantity 
receives double logarithmic corrections of infrared and colUnear origin [1], which take the form 
of double poles in dimensional regularization for the case of massless on-shell quarks studied in 
the present article. These contributions can be resummed by evolution equations in based on 
universal factorization properties of the amplitude in the relevant kinematic limit, resulting in the 
well-known exponentiation of the form factor [2-4]. So far perturbative calculations have been per- 
formed up to two loops for both the massless on-shell case [5,6] and heavy quarks [7]. Accordingly, 
the exponentiation has been studied up to the next-to-leading (ML) contributions [8-10]. 

Higher-order corrections to the quark form factor are not only of general interest in quan- 
tum field theory, but also relevant for practical applications, as this quantity contributes to phe- 
nomenologically important processes. Research in the past years has yielded dramatic progress 
in next-to-next-to-leading order (NNLO) perturbative calculations, see, for example, Ref. [11] and 
numerous references therein. This progress also led to further investigations of the general struc- 
ture of amplitudes and cross section at higher loop-orders, which in turn further stimulated the 
interest in all-order resummations. Consequently, the intimate connection between resummation 
and perturbative results at multiple loops has become much more prominent [12-14]. 

Very recently, we have presented the first complete calculation of the third-order corrections 
to a hard-scattering observable depending on a dimensionless variable, the structure function F2 in 
photon-exchange deep-inelastic scattering [15]. After exploring the consequences of that result for 
the soft-gluon threshold resummation in Ref. [16], we here present its implications on the quark 
form factor (from now on always referring to the massless on-shell case, if not explicitly indicated 
otherwise) and its resummation. We are able, after extending the two-loop form factor beyond the 
previous order s'^ in dimensional regularization, to derive the complete series of poles, e^^ . . . 
at three loops. These terms in turn provide the coefficients required to extend the exponentiation of 
the form factor to the next-to-next-to-leading (NNL) contributions which we work out explicitly. 

This article is organized as follows: In Section 2 we address the resummation of the quark form 
factor. We briefly recall the evolution equation and its solution, and present the explicit expansion 
up to four loops in terms of two perturbative functions, A{a^) (known up to three loops from the 
NNLO splitting functions [17, 18]) and G(as,e). In Section 3 we extend the two-loop form factor 
to order and extract the pole terms at three loops from our structure-function calculation [15]. 
These results are employed to extend the first- and second-order parts of the resummation function 
G to higher orders in 8, and to derive the leading-8 term at the third order in the strong coupling. 
Some first implications of these results are discussed in Section 4. Here we extend the ratio of the 
time-Uke and space-like form factors [9] to the fourth order in tts, compare to a recent result for 
9>i =A Super- Yang-Mills theory [19] and indicate appMcations on the infrared structure of massive 
gauge theories [20,21]. We briefly summarized our results in Section 5. A few technical details for 
the solution of the evolution equations in Section 2 can be found in Appendix A. Finally the break- 
up of the (8-extended) two-loop form factor into its Feynman diagrams is presented in Appendix B. 
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2. The resummation of the quark form factor 



The subject of our study are the QCD corrections to the y*^^ (or y*^^) vertex, where y* denotes a 
space-like (or time-like) photon with virtuality Q^, and ^/^ a massless external quark/ antiquark. 
Until Section 4 we will focus on the space-Uke case, thus the relevant amplitude is 

T^ = \e^(u-{^u)7(a,,Q^), (2.1) 

where the scalar function f on the right-hand side is the space-like quark form factor. This quantity 
can be calculated order by order in the strong coupling constant tts and, as mentioned above, is so 
far known to two loops [5,6]. 7 is gauge invariant, but divergent. As usual we work in dimensional 
regularization with Z) = 4 — 28, thus these divergences show up as poles in the present article. 

The exponentiation of the form factor, which extends beyond the resummation of renormal- 
ization group logarithms, is achieved by solving the well-known evolution equations [3,4,8-10] 



e'^ln!r(^as,^,ej = -K{a,,€)^ -G\^,a^,z^ . (2.2) 

Here n represents the renormalization scale, and the functions G and K are subject to the renormal- 
ization group equations [3] 

(/^^ + P(as,e) J^)G(^,as,8) = A(as), (2.3) 

(^^ + P(«s,e)3|^)^(as,e) = -A(as). (2.4) 

All infrared singularities are collected by the scale-independent function K, which in the MS 
scheme consists of a series of poles in e. The function G, on the other hand, is finite for e — > 
and includes all dependence on the scale 2^. The renormaUzation properties of G and K are both 
governed by the same anomalous dimension A, because the sum of G and ^ is a renormalization 
group invariant. This quantity is given by a power expansion in the strong couphng, for which we 
use the convention (also employed for all other expansions in tts throughout this article) 

^(«s) = l(S)'A,^f;a^A,. (2.5) 
i=\ ^^^^ i=\ 

In fact, the anomalous dimension A also occurs in many other circumstances, for instance as the 
coefficient of the 1/(1 — contribution to the AltarelU-Parisi quark-quark splitting function and 
as the anomalous dimension of a Wilson line with a cusp [22]. 

As akeady indicated by the argument of the beta function, the solution of Eqs. (2.3) and (2.4) 
requires the running coupling in D dimensions. Following Refs. [10,23] we define a(X, , e), where 
X. is a dimensionless ratio of scales like A, = Q^/i^. The resummation of the NNL contributions to 
the form factor requires the scale dependence of d to NNLO accuracy [24, 25] (see the discussion 
at the end of this section), obtained by solving 

X,^a(X,,as,e) = -ea(?L,as,8) - Poa^(X,«s,e) - Pi«^(X,«s,e) - P2a^(A.,as,e) (2.6) 
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with the boundary condition a(l,as,e) = a^. Extending the result of Ref. [23] by one order, this 
solution is given by 

.(Vs,e) = ^il-s|i^|-4|^lnX + 7)U4m^ln^Xri+7 + lr 
^ ' " ^ X\ X j X^X^o^ 'j X^X^ll \ 4 

+ I InZ Q(3 + 7)(1 -X) - 1 - 7 - ly^^ I + O {4) , (2.7) 

where we have used the abbreviations 

X = l-a,^iX-^-l), Y = '-^. (2.8) 

With Eq. (2.7) for the running coupUng, Eq. (2.3) can now be solved to the required accuracy, 

^,as,e\ = 0(1,0 (^,as,e\ ,e]+ J ^ A{a{X,as,e)) . (2.9) 

The perturbative expansion of the boundary condition G{l,d,e) can be derived by comparison to 
the fixed-order results for the form factor. 

After recursively determining (see, e.g., Ref. [10] for details) the scale-independent counter- 
term function K from Eq. (2.4), the resummed quark form factor reads 

hij (^as,^,ej = ^ I ^^*:(as,8) + G(l,a(^,a„8),s) + |^A(a(X,as,8))^ (2.10) 

with the boundary condition y (aj,0,8) = 1 [9]. After expanding the £)-dimensional couphng 
according to Eq. (2.7), In^ exhibits double logarithms of ji^ and double poles in 8, which are 
generated by the two integrations. In addition the integral over the anomalous dimension A leads 
to terms which are independent of the outer integration variable ^. These logarithmic singularities 
at ^ = are canceled by the function K order by order in the perturbative expansion. 

The well-known relation (2.10) can be employed either for a direct evaluation of the form 
factor due to the analyticity in D dimensions [10] or, by means of finite-order expansions and 
matching, for predictions of perturbative results at higher orders. Here we will focus on the latter 
issue. In particular, we will derive explicit results at three and four loops. This is done by perform- 
ing the integrations in Eq. (2.10) after inserting the perturbative expansions of all quantities. The 
resulting integrals can be evaluated using algorithms for the evaluations of nested sums [26, 27]. 
Some technical details for this step are given in Appendix A, where Eqs. (A.4)-(A.7) represent 
sample types of relevant integrals. Further details may also be found in Ref. [10]. 

It is convenient to express the loop-expanded form factor in terms of the bare (unrenormalized) 
coupling ttg instead of the renormahzed coupling tts as in Eq. (2.10). The couphngs and ttg are 
related by 

as = 2a. as, (2.11) 
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with the renormahzation constant in the MS scheme given by 

and also the bare expansion parameter normalized as = a^/ (4ji) . The perturbative expansion of 
the bare (unrenormalized) quark form factor then reads 

^F\a'„Q') = l + £(a^)' (^) '\i . (2.13) 

1=1 \M y 

In terms of the i-th order parameters A, in Eq. (2.5) and the corresponding functions G;(e), the 
expansion coefficients up to four loops read 
11 11 

2i^^^-2s' 



iFi = --^Ai---Gi, (2.14) 



^2 = ^^A? + ilAi(2Gi-|3o) + ^^(G?-A2-2|3oGi)-ilG2, (2.15) 
= -^^^1- i^^^i(^i-Po)- ii4^^i(9G?-9A2-27PoGi + 8P^) 



1 1 
T44& 

1 ^ (9Gi G2 - 4A3 - 6pi Gi - 24poG2) - ^ ^ 
1 1 4 1 1 „ . 11 



3 (3G? -9A2G1 - 18AiG2 + 4|3iAi - 18|3oG? + 16|3oA2 + 24p§Gi) 
+ ;^ (9G1G2 - 4A3 - 6P1G1 - 24P0G2) (2. 16) 



1 1 

576^5 



+ ^ ^Ai{6Gl - I8A2G1 - I8A1G2 + 8P1A1 - 45poG? + 4IP0A2 + 82p2Gi - 18p^) 
'6 

1 1 



+ ^(3^1 - I8A2G? + 9Al - 72A1G1G2 + 32A1A3 + 64piAiGi - 36poG? 

+ IOOP0A2G1 +228P0A1G2 -48poPiAi + 132p^)G? - 108pgA2 - 144p(3)Gi) 

+ ^3(-9<^?^2 + 8A3G1 +9A2G2 + 24A1G3 - 3P2A1 + 12piG? - 9P1A2 
Zoo 

+ 66P0G1G2 - 27P0A3 - 48poPiGi - 108p§G2) + {3Gl + 8G1G3 - 3A4 

96 

- 4P2G1 - I2P1G2 - 36P0G3) - ^-G4 . (2.17) 

o £ 

The three- and four-loop relations (2.16) and (2.17) are new results of the present article. Recall 
that Eqs. (2.14) - (2.17) directly refer to the bare form factor. The corresponding renormalized 
results can be derived with the help of Eqs. (2.1 1) and (2. 12). 

The £*' term of Gi, together with Po and the lowest-order anomalous dimension Ai, specify 
the two most singular terms £^2" and £~^"+^ to all orders a". Likewise, if (besides two more 
contributions to GO also the leading term of G2 and the NLO quantities Pi and A2 are known, the 
resummation fixes the first four leading poles at each order. This has been the status up to now, 
referred to as the next-to-leading (NL) contributions in Section 1. In the next section, we will 
present the leading term of G3 and the corresponding higher coefficients in the £-expansions of G2 
and G3. Together with P2 (as indicated before Eq. (2.6)) and our recent result for A3 [17], these 
results provide the NNL terms at all orders, especially fixing the e^^ and £^^ poles in Eq. (2.17). 
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3. Fixed-order results and resummation coefficients 



We now turn to the extraction of the quark form factor up to order from our third-order com- 
putation of the deep-inelastic structure functions [15]. As also discussed in Refs. [17, 18], the 
calculation has been performed via forward Compton amplitudes and the optical theorem. The cuts 
of the corresponding diagrams always include real-emission contributions, thus the purely virtual 
form-factor part caimot be directly read off at this level. Nevertheless we can reconstruct the form 
factor from our results, except (as explained below) at the highest power of £ which was consis- 
tently kept in the calculations. Consequently, we can derive all 1/e pole terms at order , since 
the forward Compton amplitudes have been computed to order 8^ for Ref. [15]. 

Our starting point for the determination of the form factor is the unrenormalized (and unfac- 
torized) partonic structure function for y*^ — >Z in the limit I, where x denotes the partonic 
Bjorken variable. Using the end-point properties of the harmonic polylogarithms [28] in which 
these results are expressed, we remove all regular contributions and only retain the singular pieces 
proportional to 6(1 —x) and the -i-distributions at order a", 

"ln*^(l-x) 



1 



k = l,...2n-l. (3.1) 

+ 



The resulting expressions are then compared to the general structure of the n-th order contribution 
F^ in terms of the /-loop form factors iF/ and the corresponding pure real-emission parts Si, 

Fo" = 6(1 -x) 

Fi" = 2^i5(l-x) + 5i 

F^ = 2^2 6(l-x) + (Ji)2 5(l-x) + 2:ri5i+52 

F^ = 2:r3 5(1 - x) + 2:Fi :F2 5(1 -x)+2iF25i+ 2^^152 + 53. (3.2) 
The x-dependence of the factors Sk is given by the Z)-dimensional -i-distributions fke defined by 

/te(x) = [(l-x)-'-^^]+ = -l5(l-x) + X^^2)/- (3.3) 

The a" contributions !Fn and Sn in Eq. (3.2) exhibit poles in e up to order e~^". The correspond- 
ing bare structure function on the other hand, only include terms up to £~", as the higher 
divergences on the right-hand sides cancel for these inclusive quantities due to the Kinoshita-Lee- 
Nauenberg theorem [29, 30]. In fact, the complete cancellation akeady occurs at the level of the 
individual diagrams of Ref. [15] for the forward Compton amplitude. 

Once the products of lower-order quantities in Eq. (3.2) have been subtracted from F^, the 
contribution of the n-loop form factor iF„ can be extracted by performing the substitution 

®o--6(l-x)-£i^®, (3.4) 

which eliminates, besides the -i-distributions, the remaining 5(1 — x) originating in the factor /„e of 
Eq. (3.3) in the purely real part Sn- However, as 5(1 —x) enters /„e with a factor l/e, this extraction 
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does not work at the highest order of £ kept in the calculation of F^. Hence, as stated above, the 
determination of the n-loop form factor iF„ to order in this approach requires the calculation of 
the bare partonic structure function to order e^"'"^ 

In addition, the subtraction of the lower-order contributions and Si with / < « in Eq. (3.2) 
requires the extension of these quantities to higher orders in 8. Specifically, the first- and second- 
order quantities are required to order 8^ and 8^ respectively, for the extraction of the pole terms 
of the three-loop form factor. These functions have been determined from the calculation of Ff 
to order 8^ and F]' to order 8^. In fact, anticipating a future extension to the finite parts of the 
three-loop form factor ^3, we have extended these calculations to one more power of 8, making use 
of the fact that the one- and two-loop integrals for the calculation of the structure functions were 
evaluated to order 8^ and 8^ anyway, see Table 3 of Ref. [15]. As a check of these new two-loop 
results (the one-loop quantities are known to all orders in 8 anyway), a separate calculation of ^2 
has been performed to order 8^ in the approach of Refs. [5, 6]. The results for the corresponding 
seven diagrams are listed in Appendix B. 

To the accuracy in 8 just discussed, the unrenormaUzed quark form factor reads, up to three 
loops in the notation of Eq. (2.13), 



2l-3i-8 + C2 + 8(-16 + ^^2 + y^3j+82(-32 + 4C2 + 7C3 + |a) 
+ 8^ (-64 + 8C2 + y C3 + '-^Q - l^2^3 + f + 8^ (- 128 + 16^2 + ^C3 

47 7 93 949 3 49 ,\\ 

+ y ^2 - 2^2^3 + y ^5 + ^^2 - y ^3 J I , (3.5) 

^2f.l .1 1/41 ^„ \ 1 /221 64„ \ 1151 17„ 
^2 = C||2^ + 6^ + ^(y-2C2j+-(^-yC3j+^ + yC2-58C3 

,0^2 /^5741 213^ 839^ 171 2 112^^ 184^ \ 2/^27911 
- 13^1 + 8 + -^2 - ^^3 - -^1 + -C2C3 - -^3) + 8^ 

I839„ 6989^ 3401„2 ^ 462^ 223^, 2608^2 

+ ^C2 - ^^3 - y^Cl + 54C2C3 - ^^5 + -^1 + -^Q 

^ ^ r 11 1 1 / 83 „ \ 1/ 4129 11„ \ 89173 



6 83 82 V 9 ^7 8 V 108 6 ^ ^7 648 
119. 467. 44„2 / 1775893 6505„ 6586„ 1891„2 

- -9-^^ + + y ^2 + £ ^-^888 lo8"^' + + ~66~^' 

89.. <. ,,<.\ 2f 33912061 146197^ 159949 „ 2639., 

- y C2C3 + 51^5 j + 8^ (-^^ - -^^2 + ^C3 + y^Cl 

397^^ 349K 809^, 569^2^ ^fll 14 1 1/^353 



U\ 7541 14 ^ 26^ / 150125 353^ 364^ 41 



3 '7 324 9 ^ 9 ' V 1944 54 ' 27 ' 30 



2/2877653 26. „ 7541. 4589. 287„2 242„\1 
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77^ 213^2\ 1/ 9073 467^ 2119^ 1461^2 214^ „ 

644^ M fll 1 1/431 \ 1 /6415 7^ 

1 /79277 1487^ 83^, \ 1 /1773839 38623^ 6703^ 

9839^2 215^^ ,,^yW ^ ^li 242 1 1 / 6521 88^ 



-72 ■ 3 ^-^-^ 'V J 81 s4 s3 V 243 27 

1 / 40289 553 1672 _ 88 2A 1 / 1870564 68497 
i2 V 243~ ~ "ST^^ ^7~^^ ~ 45 ^ e \ 2187 486~^^ 

12106„ 802„2 88„ „ 136„\) ^2 f 2 1 37 1 1 / 545 

1„ \ If 6499 133 „ 146 „ \ 1 / 138865 2849 „ 2557 

~ 3 ^^\^~lS^'^^^V^e\ ^86 5^^' + 

337^2M ^ ^ f88 1 1 /2254 16„ \ 1 /13679 316^ 

+ J / + "/"^^"^n 81 + I W - 27 ^ ^ 

256 „ \ 1 /623987 11027„ 6436„ 44„2M 2^ f 8 1 
- — ^3 ) + - ^^Y^:^ + ^^2 - -g^Ca - y ^2 j I + 8Yi4 



27 ' / £ 



188 1 1 / 124 4„ \ 1 / 49900 94„ 136„ ., 

+ :3{-^-i^^2]+-{-^T^-^^2 + ^^3] }. (3.7) 



243 V 27 9" J e\ 2187 27" 81 

Here n^r stands for the number of effectively massless quark flavours, Cf and Ca are the usual QCD 
colour factors, Cf = 4/3 and Ca = 3, and the values of Riemann's zeta function are denoted by 

Eq. (3.7) and the £^ and £^ parts of Eq. (3.6) are new results of this article. The four highest l/e 
poles of the three-loop form factor ^3 provide the first complete verification of the resummation 
of the next-to-leading contributions. With the anomalous dimensions (2.5) known up to A3, the 
remaining two poles are sufficient to fix the NNL contributions to the function G in Eq. (2.9). 
Especially, we can derive the first (e = 0) term of the third-order function Gj{e). 

Before we turn to these results we recall, for completeness, the known coefficients of the cusp 
anomalous dimension A (as). The results forAi andA2, 

Ai=4Q, A2 = 8CfQ(^^-C2)+8Cfnj(^-0 , (3.8) 

have been known for a long time [31]. The recently completed expression for A3 reads [17] 

A, = 16C,Ci - I + » C3 + ^ «) + 16C|„, (-g + 2 C3) 

+ l6C,C,nf (-^ + ^i;,-Z|;3) + i6C,„; (-1) . (3.9) 

See Refs. [32-34] for previous partial results on the «y -contributions. Very recently the C| term in 
Eq. (3.9) has been confirmed in Ref. [19], see the discussion at the end of Section 4. 
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Inserting Eqs. (3.8) and (3.9) into the resummation relations (2.14) - (2.16) and comparing 
to the explicit results (3.5) - (3.7), we obtain the following perturbative expansion of Eq. (2.9) at 

Gi = eCp + sCf (16 - 2C2) + t^Cp [32 - 3C2 - y + e^Cf ^64 - 8C2 - 14^3 

+s^C,(256-32^2-f C3-^^l + 7^2^3-f + f , O-IO) 



, 2545 44 \ / 418 8 

G2 = C|(3-24C2+48^3)+QQ ( yr + T^^"^^^V V 3^' 

^2/1 176„2\ ^^^70165 575„ 520„ 

+sC| - 1 16C2 + 120C3 + — ^1 J + Q (^^^ + — ^2 - — ^3 

176^2\ ^ ( 5813 74^ 16^ \ 2^2/^ 109 ^^^^ 
-^Cl) +eCf«/ ( - + T^O +e'c2 -437^2 + 736^3 

432^, \ , / 1547797 7297, 24958, 



2 ■ 972 27 ' 27 

653„2 356 „ „ \ 2 ^ 129389 850„ 1204„ 7„2 

+e^C| (- 1|I - ^C2 + 3614^3 + 508a - IO4C2C3 + 72^5 - 

V 5832 162 ^ 81 ^ 180 



478^ „ 1618 „3 2276 ^2\ 3 ^ / 

+— ^2C3 - 840^5 + + ) +^\Cp {- 



2628821 8405 
2916 ST^^ 



16340^ 1873^2 44^^ \ 
+^^C3 + -^a + yC2C3 + 48C5j, (3.11) 

^r,2-64r.r.-48nr.Vr2r.r^ ^096. 



G3 = C3 (^29 + 36^2 + 136^3 + ^^2 - 64^2^3 - 480^5 J + C^Q [-^ —^2 

3008„ 8„2 . ^,^.\ ^^2/^1045955 34732^ 34928^ 

+^^3 - 3 + 32C2C3 + 240^5 j + CpCl ( — ^29~ + ~8]~^' " ^7~ 



^2 176 „ „ 3826 296 „ 1232 „ 208 „ 

+ — C2C3 + 272^5 j + nfl j^- ^ + _^2 - -^C3 + ^C: 



188 

309838 11728, 1448 



^nf,c, ^-r:^ - 1^^2 + ^^3 + 

2^ / 19676 304^ 32^ \ 
+«/C^(^ + ^^^ + 27^3j. (3.12) 

As discussed at the end of Section 2, these results are sufficient to fix the next-to-next-to-leading 
contributions, i.e., the six highest poles in e, to all orders in the strong coupling. In fact, in view 
of a future extension of G3 to order £, the first- and second-order results (3.10) and (3.11) already 
transcend this accuracy by one power in e. 
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We close this section by a brief discussion of our three-loop results (3.7) and (3.12). The 
former result for the 1 /e poles of the quark form factor in massless QCD is not directly applicable to 
any physical process. For use in cross section calculations such as e^e~ ^ 2 jets at the next-to-next- 
to-next-to-leading order (N^LO), one would need the finite contribution to as well. However, 
the resulting leading term (3.12) of G3 is of immediate interest for predictions of the pole structure 
of QCD amplitudes at higher orders [13, 14] generalizing Catani's NNLO formula [12]. For the 
four-quark amplitude at A^^LO, qq qq, for instance, an explicit prediction has been derived in 
Ref. [13], for which Eqs. (3.7) and (3.12) now provide the last missing piece of information. 



4. The time-like case and non-QCD applications 

So far, our discussion has been restricted to space-like photon momenta, = —Q^ < . The 
modifications for the time-like case > are obtained by analytic continuation. For the resummed 
quark form factor in Eq. (2.10) this continuation has been discussed in Ref. [9], while the finite- 
order expansions (2.13) are transferred to > according to [5] 

-q^y'' _ (q'-y'^l n^-izw+iz) , me , ^^^^ 



/|2 J \n^) Vr(i-2/e)r(i+2/8) r(i-/8)r(i+/8) 

Of particular interest is the absolute ratio \!F {q^) / 1 {—q^) \ of the renormalized time-like and 
space-like form factors. This quantity is infrared finite and directly enters the cross section for 
Drell-Yan lepton pair production in hadronic collisions. Transforming Eqs. (2.14) - (2.17) back to 
the renormalized quantities using Eqs. (2.11) and (2.12), and then employing the analytic continu- 
ation (4.1) we obtain the expansion 

2 



7{q 



2\ 



7 i-q^) 



= 1 + as{3C2Ai} + ^Cl^f + 3C2(|3oGi + A2) | 



+ «s I \^2A\ + 3C|Ai (3poGi - Po + 3A2) + 3^2 (A3 + PiGi + 2P0G2) 

+ f C2X + ^C|A?(3PoGi - 2pg + 3A2) + ^^2'(-6P^A2 + 3^lG\ 
+ 3A^ + I2P0A1G2 + 6A1A3 + 6P1A1G1 - 5poPiAi + 6M2G1 
-6p^Gi) + 3C2(A4 + p2Gi+3|3oG3+2|3iG2)| + 0(af) (4.2) 

in terms of the couplings as(^^) = Os(~"^^) = fl^s- Note that, since this ratio is infrared finite, only 
the 8 = parts of the coefficients G, enter Eq. (4.2). Consequently, all terms contributing at the 
fourth order are now known, with the exception of the four-loop cusp anomalous dimension A4 of 
which only the small nj contribution has been derived so far [32]. 

The effect of A4 is expected to be small, therefore we can nevertheless evaluate the ratio (4.2) 
also numerically up to the fourth order, employing the [1/1] Fade estimate of Ref. [16], 

Aq,4 ^ 7849, 4313, 1553 for = 3, 4, 5 , (4.3) 
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to which we assign a conservative 50% uncertainty. Switching back to the strong coupUng ttg = 
4jias at the scale as the expansion parameter, Eq. (4.2) for n^ = A yields the numerical expansion 

2 

= l+2.094as + 5.613 a,^ + 15.70as^ + (48.63 ±0.43) (4.4) 

This result does not look like a nicely converging expansion, but so far does not exhibit a clear 
factorial growth of the higher-order coefficients either. As already pointed out in Ref. [9], the only 
genuine /-loop contribution at order al is given by the anomalous dimension A/, which in Eq. (4.4) 
contributes 24%, 7% and (2zb 1)% of the total coefficient at the second, third and fourth order, 
respectively. On the other hand, the contributions of the quantities G/_i at order are large, 
amounting to 37%, 41%, 50% at 1 = 2, 3, 4. Consequently, the higher-order (/ > 5) terms in 
Eq. (4.4) cannot be predicted quantitatively at this point. 

Exponentiations like Eq. (2.10) for the form factor iT have also been studied for electroweak 
interactions [20], where a fermion or gauge-boson mass m acts as a regulator for coUinear or in- 
frared singularities. Of course, both the counter-term function K in Eq. (2.4) and the lower integra- 
tion limit in Eq. (2. 10) are modified in this case, as they depend on the infrared sector of the theory. 
However, the leading (s = 0) term of the function G in Eq. (2.3) is independent of the regulator at 
each order in the coupling constant. This contribution entirely originates in the so-called hard re- 
gion in an expansion of the loop integrals in different regions [35]. In this region all loop momenta 
are of order Q, effectively leading to the massless case considered in Eqs. (2.16) and (3.7). 

This 'universality' implies, for instance, that Eq. (3.12) provides a prediction for the coefficient 
of ln(2^/m^) in the three-loop quantity JF3 for an Abelian gauge theory with fermion masses like 
Quantum Electrodynamics (QED) [36] after the usual identification of the colour factors. For QED, 
e.g., one has Cp = I, Ca = and Tf = 1 instead of our QCD convention Tjny = «j/2. 

Another interesting implication of Eq. (3.7) arises for maximally supersymmetric Yang-Mills 
theory (MSYM), i.e., Yang-Mills theory with 0\C = 4 supersymmetry in four dimensions. QCD 
results may be carried over to this theory using the inspired observation [37] that the MSYM results 
can be obtained from the contributions of leading transcendentality in QCD. This procedure has 
been applied to the QCD results for the three-loop anomalous dimensions of spin of leading- 
twist operators [17, 18], which were employed to extract corresponding quantities in MSYM [37]. 
Strikingly enough, the resulting MSYM anomalous dimensions completely agree with predictions 
based on integrabiUty for the planar three-loop contribution to the dilatation operator [38]. This 
agreement has been checked up to spin = 8 in Ref. [39] and is now established up to A'^ = 70 [40] 
(for a review see also Ref. [41]). 

Although no formal proof exists for the procedure of Ref. [37], it has recently been used 
in reverse, namely to predict terms of highest transcendentality in the QCD form factor. Based 
on studies of planar amplitudes in MSYM at three loops [19], where an interesting pattern of 
iteration for the four-point ampUtude has been found, both the coefficients Aily^^^M ^"^^ "^^e leading 
contribution to Gi\^^ym ^^^^ been determined for / < 3. Our new result for the three-loop form 
factor iTs in Eq. (3.7) and for coefficient G3 in Eq. (3.12) puts us in a position to check this part of 
Ref. [19] and thereby provide further evidence on the procedure of Ref. [37]. 



-10- 



The only transcendental numbers entering the results for the form factor are the values C,„ of 
Riemann's zeta function. Hence the procedure of Ref. [37] implies that, as each order in ttg, one 
keeps only the highest terms ^„ and ^; C,j with i + j = n. After the S YM identification Ca = Cf = ric 
(terms with «y do not contribute at the highest transcendentality), Eqs. (3.8) and (3.9) lead to 

176 



MSYM 



= 4nr 



A2 



MSYM 



= -Kin, 



Correspondingly, Eqs. (3.10)-(3.12) result in 



MSYM 



0, G2 



MSYM 



c ' 



G3 



A3 



MSYM 



2„3 
c ■ 



80 



MSYM 



C2C3n? + 32C5n; 



(4.5) 



(4.6) 



Both relations agree with the results of Ref. [19], and hence with the prescription of Ref. [37]. 



5. Summary 

We have derived new higher-order QCD results for the electromagnetic form factor of on-shell 
massless quarks. Specifically, we have extracted all third-order l/e pole terms in dimensional 
regularization from our recent computation of the three-loop coefficient functions for inclusive 
deep-inelastic scattering [15], supplemented by a higher-8 extension of the two-loop contributions. 
These results, together with our extension of the resummation of the form factor to the next-to- 
next-to-leading contributions, fix the six highest l/e poles to all orders. As an example, we have 
provided the explicit expression for the coefficients of . . . at four loops. 

While the pole terms of the form factor alone are not sufficient for use in other three-loop 
calculations like e~^e~ 2 jets, they do have immediate theoretical applications both for the in- 
frared structure of higher-order QCD amplitudes and for other gauge theories such as QED and 
0\C = 4 Super- Yang-Mills theory, where our results confirm a recent corresponding calculation in 
Ref. [19]. Moreover, our present results are sufficient (up to a numerically irrelevant uncertainty 
due to the unknown four-loop cusp anomalous dimension) for extending the finite absolute ratio 
of the time-like and space-like form factors, which directly enters the description of the Drell-Yan 
process, to the fourth order in ttg. 

We close by noting that the computation of the finite part of the three-loop quark form factor 
!fs by an extension of the techniques employed in this article is feasible. 
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Appendix A 

Here we give some details for the determination of the loop-expanded form factor from Eq. (2.10). 
Useful auxiliary relations are 

(A.l) 



1 


y r{n-e + i) x' 






ln*^(l-x) 




{l-x)"-^ 





(A.2) 

where Eq. (A. 1) holds for |x| < 1. The expansion of the Gamma function in powers of e for positive 
integers n reads 

nirV+e'^ = l+85i(n) + s2(5i,i(n)-52(n)) + 8^(51,1,1 (n)-5i,2(n)) 
- 52,1 (n) + 53 (n) + s'^(5i,i,i,i (n) - 5i,i,2(n) - 5i,2,i (n) 
+ 5i,3 (n) - 52,1,1 («) + S2.2 (n) + ^3,1 (n) - S^n)) + O (s^) , (A.3) 

with Smi,...,mk{^) denoting the harmonic sums [26]. Finally integrals of the following types occur: 



/ ^X-"^ = (A.4) 
J K HE 



(A.6) 



/fr«(x--ir = -i-l^(x-e-.) 

and so on, where n > 0. The integration over X and ^ in Eq. (2.10) leads to double sums which 
are readily evaluated to any finite order in a^. Also all-order analytical results for the exponent, 
cf. Ref. [10], can be obtained along these lines by employing the algorithms for the evaluation of 
nested sums [27]. 



Appendix B 

Finally we present the results for the individual Feynman diagrams, displayed in Fig 1, which 
contribute to the two-loop form factor in the approach (and notation) of Refs. [5, 6]. The diagrams 
add up to the bare quark form factor in Eq. (3.6) according to 

!r2 = 2S + QL + GL + 1QV + 1GV + C + L. (B.l) 
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Note that both the normalization in Eq. (2.13), where we have pulled out the factor [Q^ jfj?') ^, 
and our convention for 8 are different from those in Ref. [6]. 





QV CV C L 

Figure 1: The Feynman diagrams contributing to the two-loop quark form factor in the notation of Ref. [6]. 



The results for the individual diagrams to order are given by 



r ^2f 1 7 1 1/53 ^\ 355 7^ 32^ /2281 53^ 112^ 57^2 



2^14299 355„ 424„ 399„2 32„ „ 272^ 



(B.2) 



ril 14 1 1 /353 K\ 7541 14 ^ 26^ /150125 353^ 



27 30 V 11664 324 81 45 9 ^^^^ 15 ^ 

5 1 38 1 1 / 1969 5„ \ 43165 38.. 65 „ / 873877 



,/2877653 7541 4589 „ 287 26, 



242, 



GL = CpC^< - 



(B.3) 



6s3 9 82 



108 6 



648 



3888 



1969 988 41 2\ if 16929277 43165 25597 779 2 
"108"^'+ 27^^+12 ^V"^" r^3328 648~^' + ^6^^^ + 45"^^ 

65„ „ 12K M 
+ yC2C3 + — ^5 j 



(B.4) 



1 1/^11 A 1/^ 109 ^„ \ 911 91„ 59„ 
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, 6957 689 „ 296 „ 129 „2 58 „ „ ^„ \ 49639 



5'^ V 16 4 ' 3 ' 10 3 ' ' ■ '-y ■ \^ 32 

4843 1307 1267 2 293 407 281 3 M 

^^[11 1/5 1„\ 1/73 1„ 8„ \ 663 15„ 1„ 57„2 

GV=C,Q|-^ + ^^----C.j+-(-^ + 5C2-3C3j-^ + ^C2 + 5^3-^^| 



32 8 " 6 " 5"^ 3" " 5 "7 V 64 16 
1297„ 73„2 29„ „ 3„ 1343„3 128„2M /t, 

+ ^C3 + y CI - -^C2C3 + - + -^Q) I , (B.6) 

C = C^(C^-^)|i + ^ + i(16-7C2) + ^(58-16C2-^;3) +204-58^2-^^3 

53^2 (rr.. .a.y 1646^ 402^2 326^^ 842^ \ 2^4631 
- -a + e (697 - 18 IC2 - ^C3 - ^Cl + ^^2^3 - j + 

1141 6293 1744 2,836 2708 1399 3 4274 2\ 1 

" " " + —^2^3 - + ^^2 + ^^3 J I , (B.7) 

r ^2(1 2 1/17 ^\ 1/101 46^ \ 631 35^ 152^ 103^2 

/3941 335^ 439^ 159^2 98^^ 598 ^ \ ,/24495 2573^ 
+ e (^ - — ^2 + — ^3 + — ^1 - y ^2^3 + — ^5 j + e- (^ - ^^2 

2065„ 1839„2 152 „ „ 1976„ 2847„3 1318„2M o^ 

+ + - -^^^^ + + " / • ^^-'^ 

The loop integrations have been reduced with integration-by-parts identities [44, 45] to so-called 
master integrals. This step has been automatized in Ref. [46]. All master integrals except the 
basic non-planar triangle can be expressed in terms of Gamma functions, thus they can be readily 
expanded to any order in s. The non-planar triangle (see, e.g., Ref. [47]) can be written as a double 
sum over Gamma functions. After expansion, the sums can be solved in terms of the Riemann zeta 
function to any order in e using the algorithms for harmonic sums [26] coded, as all our symbolic 
manipulations, in FORM [48]. 
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